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Abstract: Infiltration is an important part of the hydrological cycle, and it is one of the main abstractions 
accounted for in the rainfall-runoff modeling. The main purpose of this study is to compare the 
infiltration models that were used to assess the infiltration rate of the Mitidja Plain in Algeria. Field 
infiltration tests were conducted at 40 different sites using a double ring infiltrometer. Five statistical 
comparison criteria including root mean squared error (RMSE), normalized root mean squared error 
(МКМР), coefficient of correlation (CC), Nash-Sutcliffe efficiency (NSE), and Kling-Gupta efficiency 
(КСЕ) were used to determine the best performing infiltration model and to confirm anomalies between 
predicted and observed values. Then we evaluated performance of five models (1.е., the Philip model, 
Kostiakov model, Modified Kostiakov model, Novel model, and Horton model) in simulating the 
infiltration process based on the adjusted performance parameters cited above. Results indicated that the 
Novel model had the best simulated water infiltration process in the Mitidja Plain in Algeria. However, the 
Philip model was the weakest to simulate the infiltration process. The conclusion of this study can be 
useful for estimating infiltration rate at various sites using a Novel model when measured infiltration data 
are not available and are useful for planning and managing water resources in the study area. 
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1 Introduction 


Infiltration refers to the process of water entering the soil or rock layers from the surface. This 
phenomenon can occur naturally, for example, through precipitation, or it can be facilitated through 
human activities such as irrigation or watering plants (Angelaki et al., 2021). Infiltration of water 
into the unsaturated zone is one of the most important components of the hydrologic cycle (Iovino 
et al., 2017). Over the last few decades, several physical and empirical models have been developed 
to quantify the infiltration process, particularly for one-dimensional vertical infiltration (Mishra et 
al, 2003). The Richards (1931) equation is one of the classic physical models for calculating 
semi-infinite, isotropic, and stable soil infiltration processes. Infiltration rate depends on the nature 
of the soil (proportions of sand, silt, and clay) and the mineralogy of clay. Water moves faster 
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through the large pores of sandy soil than through the small pores of clay, especially when the clay 
is compacted and has little or no structure or aggregation (Schumann, 1998). These soil infiltration 
properties are determined when mathematically fitting an infiltration model (Thomas et al., 2020). 
But not all infiltration models are suitable for all soils (Oku and Aiyelari, 2011). Many researchers 
compare the accuracy of different models by comparing calculated and observed infiltration rates. 
Under different conditions, a particular model will show better prediction than others. But so far, 
there is no explicit mention of which model gives the best prediction (Kumar and Sihag, 2019). 
Many studies have been done for comparison of different infiltration models. The Kostiakov and 
modified Kostiakov models were the most efficient in estimating the infiltration rate and cumulative 
infiltration, respectively. The Kostiakov model was most suitable for assessing plain in the Okhlal 
landfill, Delhi, India. The ability of new infiltration model to estimate infiltration from actual field 
data was investigated in comparison with the Philip, Kostiakov, Horton, SCS, and modified 
Kostiakov models under similar conditions (Jain and Chakma, 2023). 

Infiltration plays a crucial role in irrigation because it determines how efficiently water is 
absorbed by the soil and made available to crops. If the infiltration rate is too slow, water can pool 
on the surface and cause runoff, which can lead to soil erosion and water loss (Adeniji et al., 
2013). On the other hand, if the infiltration rate is too fast, water may not be retained in the soil 
long enough for plants to absorb it, resulting in water stress for crops. By measuring and 
managing infiltration, farmers can optimize their irrigation practices and ensure that water is used 
effectively and efficiently to promote crop growth and yield (Thomas et al., 2020). 

This research aims to experimentally estimate and compare soil infiltration rate in the Mitidja 
Plain in Algeria and five different infiltration models, i.e., Novel, Kostiakov, modified Kostiakov, 
Philip, and Horton. For this study, a number of specific sites within the region were selected and 
samples were collected from those sites. The infiltration rate was determined by a double ring 
infiltrometer. Statistical analyses were performed to compare infiltration models and infiltration 
rate observed. This study is the first step in a project aimed to improve irrigation water use and 
evaluate irrigation systems in the Mitidja Plain. 

The Mitidja Plain is a traditional irrigation area that is originally composed of colonial lands 
devoted to viticulture and arboriculture. The Mitidja Plain is considered as one of the most fertile 
areas of Algeria and of the entire world due to its agricultural production. Water consumption in 
the agricultural sector of the plain represents almost 70% of the total volume of mobilized water, 
which is a considerable volume of water at current situation of water resources (Meddi and 
Eslamian, 2020), as in other regions of the world. This situation is prompting decision-makers to 
find more economical techniques to cope with climate change and the overexploitation of 
groundwater. In order to develop innovative techniques and approaches, it is essential to study the 
evolution of infiltration into the ground in order to set up intelligent control and irrigation systems. 
Understanding of infiltration in agricultural soils is an important phase in the design, evaluation, 
and management of furrow irrigation (Esfandiari and Maheshwari, 1997). 


2 Materials and methods 


2.4 Study area 


The Mitidja Plain is located in the north of Algeria, facing the Mediterranean Sea in the northeast 
(36?22'30"—36?48'58"N, 02?32'07"—03?28'14"E; 300 m а.5.1.). It is one of the main agricultural 
plains in the Mediterranean Basin. It has a total area of nearly 1700 km? with agricultural areas 
ranging from 120x10? to 130x10? hm2. It extends from the Wadi Nador in the west, bordered by 
the Dahra Mountains, east of the Oued Boudouaou, from the foothills of the Atlas Blidéen to the 
Sahel Mountains in the south, and the Bouzaréah Massif in the north (Bouziane et al., 2021). 
Administratively, it includes four cities: Algiers, Blida, Boumerdes, and Tipaza. A Mediterranean 
climate prevails over the Mitidja Plain. Summers are hot and dry, and winters are rainy. The 
average winter temperature is 6?C, the average summer temperature is 33?C, and the average 
annual precipitation is 500—700 mm. The wet season is from October to May and the dry season is 
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from June to September, and the estimated potential evapotranspiration in western and eastern 
Mitidja is about 1326 and 1271 mm, respectively (Boufekane et al., 2020). The hydrology of the 
study area is characterized by intense channel densities of four major rivers: Wadi Legaia, Wadi 
Hamiz, Wadi El Harrach in eastern Mitidja, and Wadi Mazafran in western Mitidja. The 
tributaries of these rivers originate from the Blidean Atlas and have considerable north-south 
flow. Figure 1 shows the experimental sites and soil texture. This area is characterized by fine and 
fertile soil. The most common soil texture is clay, and soils are dominated by alluvial and 
colluvial soils. Most farms are on alluvial soil. 
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Fig. 1 Study area with selected sites and soil texture (Boufekane et al., 2021; Bouziane et al., 2021) 


2.22 Measurement of infiltration rate 


To measure the infiltration rate, we used a double ring infiltrometer. The double ring infiltrometer 
is a device used to measure the infiltration rate of water into soil. It consists of two concentric 
rings 250 mm deep, the inner ring (110 mm) being smaller than the outer ring (320 mm), with the 
space between the two rings creating a circular trench (Fig. 2). The rings are made of metal and 
are driven into the soil to 40 mm depth, leaving only the top edge of the rings exposed, both rings 
were filled with water to the equal depth, and the initial reading of the water level was noted. The 
depth of water in the inner ring of infiltrometer was noted at regular time intervals (1 min) until 
the rate of infiltration became steady. The infiltration rate is typically calculated as the volume of 
water that infiltrates into the soil per unit of time and surface area of the soil. The double ring 
infiltrometer is a useful tool for assessing the infiltration capacity of soil and can be used to 
evaluate the effectiveness of soil conservation measures, irrigation practices, and other land 
management practices (Soil Science Equipment | Double Ring Infiltrometer Surface Measure, 2020). 
It is particularly useful in agricultural area where efficient water use is important. 


2.3 Experimental data 


Soil samples were collected with core samplers to an average depth of 50 cm, and disturbed soil 
samples (40 points) were collected in the field at the site and air-dried, crushed, and passed 
through a 2-mm sieve to obtain soil fractions for soil texture determination. We used a textural 
classification triangle to classify the soil based on the results obtained from the analysis. 
Experimental tests were conducted individually for 40 sites (experimental samples; Fig. 1). The 
result of analysis of soil physical properties and infiltration rate (initial and final) of the study area 
at each site is presented in Table 1. The results showed that the texture of sampled sites was 
predominantly clayey. The initial infiltration rate varies from 11.70 to 330.00 mm/h and the range 
of final infiltration rates was from 3.71 to 138.42 mm/h. Variations in infiltration rates are favored 
by the importance of root system and animals digging into the soil, insufficient pre-wetting, and 
soil disturbance by double rings infiltrometer. 
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Fig.2 Double ring infiltrometer 
2.4 Infiltration models and parameters 


In this study, infiltration data were fitted to the infiltration models of Kostiakov, Philip, Horton, 
modified Kostiakov, and Novel. 
2.4.4  Kostiakov model 
Kostiakov model is commonly used in irrigation engineering and is derived using data observed 
in the field or in the laboratory (Parhi et al., 2007). Kostiakov (1932) proposed the following 
best-fit infiltration equation: 

fü)-axt ^, (1) 
where ft) is the infiltration rate (mm/h); a is the coefficient; ¢ is the time (h); and b is the 
exponent. Both the parameters rely on soil type, initial moisture content, rainfall rate, and 
vegetative cover, and these can be determined empirically. 
2.4.2 Modified Kostiakov model 


Modified Kostiakov equation includes the term steady infiltration capacity (fc; mm/h). The 
parameter fc is included because the infiltration rate decreases as more water penetrates the soil, 
eventually reaching a constant or final infiltration rate (Vand et al., 2018). The modified 
Kostiakov equation is expressed as follows: 

Kt)=axt + fe. (2) 
2.4.3 Philip model 
Philip and De Vries (1957) proposed a model based on a semi-analytical solution of the Richards 
flow equation (Richards, 1931). In this solution, Philip and De Vries (1957) used a time series to 
solve the Richards equation. Neglecting the higher order terms, we derived the results of a 
two-term infiltration equation in which both parameters have physical meanings (Zakwan et al., 
2016). 


уже +K, (3) 


where 5 is the Sorptivity parameter that is the function of soil matrix forces (mm/h 95); and K is 
the soil parameter related to transmission of water through soil or gravity force. 

2.4.4 Horton model 

Horton equation is one of the most popular empirical models simulating water infiltration into 
soils (Dahak et al., 2022). Infiltration equation is a three-parameter equation that is expressed as 
follows: 


f) fe o-fe)xe “, (4) 
where fo is the initial infiltration rate (mm/h); fc is the final infiltration rate (mm/h); e is the 
mathematical constant; and K is the rate constant in the dimension of time t (1/4). 
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Table 1 Details of the initial and final infiltration rates (f(t)), and soil properties of the study area 
Site Coarse Fine silt Clay Fine sand ирети Initial ДӘ) at = Final/steady 
silt (%) (99) (99) (9) sand (96) 1 min (mm/h) (t) (mm/h) 
Pl 5.0 12.5 50.0 17.5 20.0 Sandy clay 240.0 30.4 
P2 12.0 13.5 40.0 12.6 32.5 Sandy clay 180.0 48.4 
P3 5.3 22 50.0 17.0 30.3 Sandy clay 77.6 30.0 
P4 12.5 2:9 55.0 37.5 7.5 Sandy clay 60.0 13.8 
P5 12.0 15.0 40.0 12.3 26.3 Clay-limono 120.0 62.8 
P6 15.0 12.5 41.5 18.3 22.8 Clay-limono 180.0 11.7 
P7 50.0 5.0 34.5 des 7.5 Clay-limono 65.5 48.0 
P8 20.0 20.0 55.0 13.9 14.8 Clay-limono 60.0 13.7 
P9 15.0 15.0 10.0 51.0 22.6 Sandy loam 120.0 22.4 
P10 1.0 4.5 35.0 12.0 65.0 Sandy loam 330.0 138.4 
РИ 10.0 55.5 25.0 6.2 77 Fine silt 90.0 61.2 
Р12 34.0 25.0 17.0 19.3 4.8 Fine silt 48.0 12.0 
P13 3.6 70.0 5.8 13.2 11.8 Fine silt 90.0 51.0 
P14 45.0 25.0 9.4 18.1 5.7 Fine silt 180.0 29.7 
P15 15.0 56.5 15.0 11.2 5:3 Fine-clay-limono 180.0 29.6 
P16 75.0 7.3 2.8 9.7 3.7 Fine-clay-limono 105.0 60.0 
P17 35.0 35.0 20.0 10.8 7.3 Fine-clay-limono 120.0 15.5 
P18 9.0 75.0 11.5 3.7 2.5 Fine-clay-limono 30.0 6.0 
P19 5.0 15.0 42.5 16.2 322 Sandy clay-loam 180.0 15.0 
P20 10.0 5.0 52.4 19.3 16.7 Sandy clay-loam 24.0 4.6 
P21 7.0 10.0 41.5 15.2 26.8 Sandy clay-loam 36.0 4.4 
Р22 4.8 2.7 40.0 19.2 40.0 Sandy clay-loam 120.0 43.3 
P23 5.0 7.5 75.0 8.4 12.1 Clay 180.0 38.0 
P24 0.5 15.0 75.0 5.1 3.8 Clay 168.0 93.8 
P25 30.0 5.0 55.0 4.9 14.2 Clay 60.0 14.4 
P26 5.0 16.0 75.0 4.3 2.8 Clay 180.0 88.8 
P27 9.0 9.0 63.0 7.6 20.6 Clay 11.7 6.0 
P28 10.0 22.5 70.0 5.7 2.6 Clay 16.4 12.0 
P29 20.0 10.0 45.0 22.3 4.5 Clay 120.0 56.7 
P30 15.0 15.0 65.0 9.0 3.7 Fine-clay-loam 180.0 43.0 
P31 35.0 10.0 40.0 5.5 10.4 Fine-clay-loam 60.0 9.6 
P32 32.0 12.0 45.3 7.2 4.4 Fine-clay-loam 15.0 4.0 
P33 20.0 14.0 50.6 Tal 8.4 Fine-clay-loam 60.0 3.7 
P34 19.0 21.3 59.7 4.4 3.8 Clay-loam 60.0 48.0 
P35 35.0 2.8 52.1 10.5 0.9 Clay-loam 120.0 34.8 
P36 21.0 5.0 60.6 10.9 2.5 Clay-loam 60.0 5.9 
P37 45.0 45.0 4.4 7.0 0.5 Very fine silt 80.0 55.4 
P38 55.0 30.0 5.3 11.5 2.0 Very fine silt 60.0 26.7 
P39 79.0 15.6 0.5 4.8 3.5 Very fine silt 27.3 8.9 
P40 16.0 68.5 5.7 5.6 5.4 Very fine silt 60.0 20.8 
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2.4.5 Novel model 


A typical non-linear regression model developed from the modified Kostiakov model (Sihag et 
al., 2017) was adopted, as shown below: 


ДӨ=ахг?+с}, (5) 
where c 1s the coefficient. 


2.5 Performance evaluation parameters 


2.5.1 RMSE 

RMSE is one of the most commonly used measures of prediction quality. It shows how far the 
predictions deviate from actual measured values (Bennett et al., 2013). RMSE is a good measure 
of accuracy, but only for comparing the prediction errors of different models or model 
configurations for a particular variable, not between variables. 


nase [taa] (6) 


i-0 


where a; is the simulated value; 5; is the observed values of infiltration rate; and и is the number 
of observations. 

2.5.2 NSE 

NSE criterion of Nash and Sutcliffe (1970) was introduced as a measure of model performance. 
NSE criterion better reflects the desirable and undesirable properties of model of interest and 
increases as the quality of model improves. NSE can be defined as follows: 


> (O -ay) 
25 (b -hy) 


where b the average of observations values. For nonlinear models, NSE can be negative. In 
fact, it ranges from —oo to 1. As usual (Zeybek, 2018), researchers look for NSE values close to 1. 
Negative values indicate unacceptable model performance. This is one of the recommended 
model performance metrics in hydrology as well as many other applications. 


2.5.3 KGE 
KGE fills several gaps in NSE and is increasingly used for model calibration and evaluation. 

KGE -1-(r-1Y « (8-1Y +(C-1Y, (8) 
where r is the linear correlation coefficient between observed and computed values; В is the ratio 
between SD of simulated flows and SD of observed flows; and C is the ratio between mean 


simulated flow and mean observed flow (Knoben et al., 2019). This coefficient is more accurate 
than NSE. 


2.5.4 CC 


CC is a statistical measure that calculates the strength of relationship between movements of two 
variables (Lei et al., 2020). CC is calculated as follows: 


Ue У. (ав) 3 (0070) 
| Y ec (Ya | | | 
2.5.5 NRMSE 


NRMSE makes И easy to compare models at different scales. NRMSE refers to RMSE of 
observed range of variables (Tadj et al., 2014). Therefore, NRMSE can be interpreted as part of 
the total range normally resolved by the model. 


NSE=1 


(7) 


C 


(9) 


zd 
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NRMSE = | > . (10) 


2.5.6 Percent bias (PBIAS) 


PBIAS measures the average tendency of simulated values to be larger or smaller than observed 
values (Golmohammadi et al., 2014). 


1 
PBIAS=100— ;-bh). П 
> (a-h) (11) 


i=0 


3 Results 


3.1 Model parameters 


To estimate infiltration rates in our study area, we opted for the models of the Novel, Kostiakov, 
modified Kostiakov, Philip and Horton. Observed field data on infiltration rates were used to 
evaluate these models, and numerical values for the infiltration equations parameters were 
obtained. The obtained model parameters and statistical parameters were done using the Rtop 
package of R statistical software. 

A comparative analysis of soil water-infiltration characteristics for ten different classes of soil 
was performed. The estimated parameter values for the five selected infiltration models are listed 
in Table 2. The values of 5; for the Novel, Kostiakov, and modified Kostiakov infiltration models 
were of significant similarity across the ten soil types. In contrast, parameters f: were completely 
different in the Novel, Horton, and modified Kostiakov models due to other parameters such as c 
in the Novel model and А in the Horton model. These empirical parameters have a soil physical 
meaning. A size comparison of each of the model parameters indicated that there is no 
relationship between the parameter k in the Horton model and Philip model. The parameter a; in 
the modified Kostiakov model is similar to that in the Novel model. However, this is not the case 
in the Kostiakov model. 

Individual infiltration curves were generated from experimental tests. Results of the field 
infiltration rates and simulated infiltration rates for the five models for each soil class are 
presented in Figure 3. The graphics showed that the observed infiltration data took 35 min to 3 h 
to reach the steady infiltration rate. 

The Philip model had the highest deviations because it overestimated the data measured at 
sampling points P1 and P15. The Horton model had the second lowest performance after the 
Philip model. This may be due to the fact that their parameters do not have a consistent physical 
interpretation and that the process involved in evaluating the parameters can be very sensitive to 
approximation errors. 

For clay-sandy, fine-clay-loam, clay-loam, and clay soils, the Horton and Philip models 
deviated from the observed data, while the modified Kostiakov and Novel models provided a 
better fit. For fine-silty and very fine-silty soils, it was difficult to decide which model gave the 
best performance, so the performance parameters were used to make a better comparison. 
However, the infiltration rate estimated by the Novel model was the most effective in predicting 
the fit of the experimental data measured at all sampling points. 


3.2 Performance parameters 


Statistical criteria such as RMSE, NRMSE, CC, NSE, and KGE were used to evaluate infiltration 
models. We selected the best model based on the following criteria: the maximum CC, the 
minimum RMSE, the minimum NRMSE, the maximum NSE, and the maximum KGE. The 
results shown in the box plot (Fig. 4) were computed for different models and the values are 
shown in Table 3. 

The KGE and NSE values were close to 1.00, with intervals of 0.90—0.99 and 0.80-0.98, 
respectively, for the Horton model in sandy clay, clay-loam, and sandy loam soils. This indicates 
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Table2 Parameters of estimated infiltration models 


Philip Horton Kostiakov Modified Kosiakov Novel 


Site 
K S fo te K di bi di bi te qi b, c £ 


Pl 6.12 49.14 352.80 41.14 30.05 29.77 0.46 6.76 0.82 28.44 3.65 0.98 5.74 5.83 
P2 38.76 47.40 158.99 54.54 3.37 60.79 0.29 60.79 0.29  2.3E-03 60.79 0.29 1.7E-05 4.61 
РЗ 40.61 16.55 67.79 0.11 0.28 49.09 0.15 49.07 0.15 4.2Е-04 49.08 0.15 10.90 6.3E-06 
Р4 12.66 15.72 41.31 8.31 0.97 20.50 0.28 20.50 0.28  8.6E-05 20.50 0.28 1.8Е-05 1.36 
Р5 59.44 12.55 178.43 68.92 44.54 63.57 0.11 1.56 0.85 65.48 0.44 1.17 28.13 2.40 
P6 3.6Е-06 37.22 206.60 17.04 14.35 10.31 0.72 10.31 0.72  7.7E-05 10.31 0.72 6.9E-06 9.26 
Р7 53.64 7.16 66.29 28.01 0.30 56.95 0.07 56.96 0.07  2.1E-04 56.95 0.07 8.29 1.6Е-05 
Р8 7.09 11.63 102.15 18.21 42.37 11.48 0.36 2.03 0.74 12.83 0.52 1.08 8.19 1.95 
P9 9.31 37.69 120.40 26.07 5.14 27.64 0.41 27.65 0.41 1.6E-03 27.64 0.41 8.9E-05 1.83 
P10 125.73 66.06 304.21 100.40 2.66 142.75 0.22 142.75 0.22  6.0E-04 142.75 0.22 5.3Е-04 0.60 
Pll 61.96 11.72 9001 55.15 1.74 64.99 0.11 65.00 0.11 1.8Е-03 65.00 0.11 5.12 1.8E-05 
P12 7.56 11.02 55.70 17.11 15.58 11.31 0.36 11.31 0.36  2.6E-03 11.31 0.36 1.04 8.8E-06 
P13 53.89 12.20 77.66 50.25 1.14 59.44 0.10 5941 0.10  1.6E-03 59.42 0.10 6.36 4.9E-05 
Р14 17.03 43.38 144.99 35.62 5.07 38.39 0.36 35.14 0.38 3.29 36.76 0.37 2.76 0.60 
P15 4.03 43.71 239.41 33.61 19.55 25.92 0.46 7.73 0.76 21.43 6.66 0.80 6.02 3.82 
Р16 81.23 14.41 112.55 8.9Е-04 0.39 85.99 0.10 86.02 0.10 6.9Е-04 86.00 0.10 7.9Е-05 0.01 
P17 1.8Е-04 35.98 135.59 20.06 9.28 16.51 0.53 16.51 0.53 9.2Е-05 16.51 0.53 9.9E-06 0.97 
P18 405 983 2794 718 3.02 891 0.36 8.91 0.36 2.3Е-04 8.91 0.36 3.4E-06 1.14 
P19 7.4Е-06 47.38 190.14 20.59 8.86 17.55 0.61 17.55 0.61  3.4E-05 17.55 0.61 0.61 1.7E-05 
P20 1.71 6.16 34.99 6.87 23.05 4.45 0.42 4.45 0.43 0.36 4.45 0.42 1.56 1.4Е-05 
Р21 3.2Е-05 11.05 34.63 4.20 5.73 5.04 0.53 5.04 0.53 5.5E-05 5.04 0.53 9.7Е-07 0.20 
Р22 35.07 23.69 133.99 49.10 11.71 45.27 0.22 9.16 0.56 38.12 9.13 0.56 1.44 26.45 
Р23 29.86 35.79 162.58 47.80 8.62 47.29 0.28 18.78 0.49 28.91 16.62 0.52 9.34 3.34 
P24 98.32 31.31 170.33 60.20 1.21 107.97 0.15 107.99 0.15 3.2Е-04 107.97 0.15 0.01 9.1Е-05 
P25 12.95 8.46 77.16 1622 25.17 17.83 0.17 0.87 0.95 15.56 0.78 0.98 5.97 2.62 
Р26 91.15 25.56 159.85 96.29 3.53 98.41 0.14 98.40 0.14 1.4E-03 98.39 0.14 2.47 1.3Е-04 
Р27 8.35 1.94 1113 6.74 0.37 9.50 0.09 9.50 0.09 2.1E-05 9.50 0.09 1.4E-05 6.22 
Р28 12.75 2.51 16.65 6.51 0.26 14.14 0.08 14.14 0.08  1.1E-04 14.15 0.08 0.12 3.9Е-05 
Р29 56.29 15.76 101.44 63.08 4.77 62.45 0.13 25.38 0.25 37.77 1423 0.37 9.12 5.41 
P30 28.13 51.71 167.76 47.58 4.01 51.83 0.34 51.83 0.34  5.3E-05 51.83 0.34 7.4E-05 2.34 


P31 1.2Е-04 13.65 110.15 14.47 44.32 6.49 0.51 2.43 0.74 6.65 0.75 1.02 2.42 4.43 
P32 268 383 17.21 476 9.32 448 0.31 2.63 0.42 1.92 2.63 0.42 0.97 1.98 
P33 1.1Е-05 11.56 94.09 643 30.41 7.61 0.38 1.96 0.80 3.47 1.47 0.88 1.04 3.92 
P34 4828 4.90 61.83 28.67 1.09 48.08 0.07 48.09 0.07 3.7Е-04 48.08 0.07 3.2E-05 2.96 
P35 27.85 24.97 105.82 40.37 6.09 38.05 0.27 32.78 0.29 5.54 32.78 0.29 2.76 2.01 
P36 3.72 10.57 85.78 7.53 25.26 9.57 0.34 2.19 0.76 6.06 1.91 0.80 3.81 1.64 


Р37 46.37 48.22 175.34 64.38 3.92 67.98 0.28 67.98 0.28 8.0Е-05 67.99 0.28 5.47 6.5Е-05 
Р38 22.93 15.93 65.68 29.08 3.99 30.30 0.22 30.30 022 1.4E-04 30.30 0.22 1.1E-04 2.53 
P39 694 593 2889 1107 11.03 8.60 0.29 8.59 0.29 1.2E-04 8.60 0.29 4.97 1.7Е-05 
Р40 17.73 8.68 6942 22.8 19.95 22.4 0.16 2.56 0.64 19.62 1.04 0.88 17.76 1.18 


Note: all parameters of these models are explained in Section 2.4. 
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Fig.3 Comparison of observed infiltration rate with estimated infiltration rate from various models in different 
sites. (a), P1 site; (b), P6 site; (c), P9 site; (d), P12 site; (e), P15 site; (f), P20 site; (g), P25 site; (h), P33 site; (1), 
P36 site; (1), P39 site. 


that the Horton model can predict the infiltration rate for these soil classes. For the remaining soil 
classes, the Novel model had the highest values, between 0.85 and 0.99, for both parameters. 

The NRMSE values for the Philip model were unacceptable, ranging from 0.03 to 0.35. In 
contrast, the Novel model had the best results, with NRMSE values ranging from 0.01 to 0.19. 
The Kostiakov and modified Kostiakov models produced good results compared with the Philip 
model, with the KGE and NSE values ranging from 0.57 to 0.99 and 0.15 to 0.99, respectively, 
for both models. The NRMSE values for the Kostiakov and modified Kostiakov models ranged 
from 0.03 to 0.43 and 0.02 to 0.19, respectively. 

The KGE parameter in Figure 4a showed that the Novel and modified Kostiakov models had 
the median value closest to the ideal value of 1.00, with a range of 0.91—0.99 and five outliers for 
both models. In contrast, the Philip model had the largest value ranging from 0.84 to 0.93, with 
three outliers. The Horton model had a uniform distribution ranging from 0.91 to 0.99, with two 
outliers. 

For the NSE parameter in Figure 4b, the modified Kostiakov and Novel models gave the same 
results, with a range of 0.85—0.98 and five outliers. The Kostiakov model had the largest value 
ranging from 0.75 to 0.97, with only one outlier. Like the KGE parameter, the Horton model had 
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a uniform distribution ranging from 0.83 to 0.97, with two outliers. Finally, for the NRMSE 
parameter in Figure 4d, the Novel and modified Kostiakov models gave the best results 
compared with the other infiltration models, approaching the ideal value of 0.00 with a range of 
0.02—0.14. The Philip model had the highest uncertainty, with a range of 0.06—0.16, followed by 
the Kostiakov model (0.05—0.12) and the Horton model (0.05—0.12), in terms of estimating 
infiltration rate values. The Novel and modified Kostiakov models performed the best in terms 
of all three statistical criteria (KGE, NSE, and NRMSE) and the Philip model performed the 
worst. 

The Novel model explained almost all the variability in the data, as indicated by its high KGE 
and NSE values and low NRMSE value. A comparison of the statistical parameters shows that the 
Novel model fits well with the measured data for all types of soil except sandy clay, clay-loam, 
and sandy loam soils, where the Horton model performed better. 

Figure 5 showed the PBIAS errors for different infiltration models. The Novel and modified 
Kostiakov models had the lowest PBIAS errors, ranging from —12% to 1296, with most of the 
values concentrated around 0%. The Philip model had the largest PBIAS errors, ranging from 
—20% to 20%, with a significant number of values concentrated around 5%. The Kostiakov model 
had an asymmetric PBIAS distribution, ranging from —27% to 15%, similar to the Horton model 
(-15%—20%). Therefore, we can conclude that the Novel model has the best fit to the measured 
data, as it is the closest to the ideal value of 096 for PBIAS error. This is likely because the Novel 
model is developed from the modified Kostiakov model. The Novel model outperformed the other 
infiltration models in terms of statistical criteria and PBIAS error. It is therefore the most suitable 
model for predicting infiltration rates for a variety of soil types. This means that the Novel model 
explained almost all the variability in the data. A comparison of the statistical parameters shows that 
the Novel model fits well with the measured data. The Novel model gave the best fit to the measured 
infiltration rate for all types of soil except for sandy clay, clay-limono, and sandy loam soils. 

Figure 5 showed that the Novel and modified Kostiakov models gave better results because they 
had the least errors such that the PBIAS was between —12% and 1296, and had the most density 096, 
whereas the Philip model had a large variation ranging from —20% to 20% with an important 
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Fig. 4 Comparision of infiltration models based on KGE (a), NSE (b), CC (c), and NRMSE (d) statistics. Boxes 
indicate the IQR (interquartile range, 75'^ to 25'^ of the data). The median value is shown as a line within the box. 
Outlier is shown as circle. Lines extend to the most extreme value within 1.5xIQR. KGE, Kling-Gupta efficiency; 
NSE, Nash-Sutcliffe efficiency; CC, coefficient of correlation; NRMSE, normalized root mean squared error; 
M-Kostiakov, modified Kostiakov model. 
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Fig. 5 Distribution of PBIAS (percent bias) of different infiltration models. (a), Horton model; (b), Kostiakov 
model; (c), Philip model; (d), M-Kostiakov (modified Kostiakov) model; (e), Novel model. 


density around 5%. The Kostiakov model had asymmetric PBIAS plot ranging from —27% to 1596, 
and Horton model from —15% to 2096. So, we can conclude that the Novel model has the best fit in 
view of the fact that the Novel model is developed from the modified Kostiakov model. 


4 Discussion 


In the context of our current research, we applied infiltration models across ten distinct soil types 
with the primary aim of evaluating the performance of five infiltration models when compared 
with observed data. This comprehensive analysis allowed us to gain valuable insights into the 
effectiveness of these models in accurately predicting infiltration characteristics in a diverse range 
of soils. The results of our study, which are presented in Table 3, and Figures 4 and 5, 
demonstrated the superiority of the Novel model in replicating observed experimental data. This 
model demonstrated remarkable success in forecasting infiltration behavior in seven different soil 
types, with exceptional accuracy. It is important to note that, while the Novel model excelled for 
the majority of soil types, three exceptions were identified, namely sandy clay, clay-limono, and 
sandy loam soils, which were best characterized by the Horton model. This finding aligns with the 
results of Singh et al. (2018), who predicted cumulative infiltration rate under inceptisols in 
tropical dry brown soils and found the Novel and modified Kostiakov models to be more suitable 


Amina MAZIGHI et al.: Estimation and inter-comparison of infiltration models in the... 1487 


than the Kostiakov and SCS (soil conservation service) models. Furthermore, our assessment 
encompassed a broader scope, extending to the region of Lafia in North Central Nigeria, where 
four distinct infiltration models were evaluated, namely the Kostiakov, modified Kostiakov, 
Philip, and Horton. In this specific context, the Horton model emerged as the most suitable fit for 
measured cumulative infiltration rate, as affirmed by Ogbe et al. (2011). This particular finding 
holds significance for the practical application of infiltration modeling in that region. Additionally, 
Wang et al. (2017) concluded that the Mezencev and Horton models offer a more accurate 
description of the relationship between cumulative infiltration and infiltration time when biochar 
is introduced into the soil. Moreover, as deduced by Duan et al. (2010), the Mezencev and the 
Horton models surpassed the other three infiltration models in terms of performance. Their 
researches provided the most accurate depiction of the intricate relationship between cumulative 
infiltration and infiltration time within the lawn soils. This finding underscores the superior 
capabilities of the Mezencev and Horton models, shedding light on their ability to offer a more 
precise representation of infiltration rate. 

It should be noted that the result of Duan et al. (2010) closely mirrors our findings for sandy 
clay, clay-limono, and sandy loam soils, where the Horton model exhibited a superior fit when 
compared with the Novel model. This intriguing consistency in results further gains support from 
parallel experiments conducted at the Agricultural Research Farm of the Bahauddin Zakariya 
University in Multan, Pakistan, as detailed by Farid et al. (2019). Their research involved the 
assessment of three distinct infiltration models, namely the Kostiakov, Philip, and Horton models. 
Their results reinforce the effectiveness of the Horton model in accurately describing infiltration 
behavior. Moreover, Rasool et al. (2021) found that the modified Kostiakov model consistently 
outperforms other models in terms of infiltration predictions. Our study also confirmed its 
modeling power except for the Horton model. These results underscore the necessity of 
thoroughly assessing and documenting the performance of various infiltration models at specific 
locations. Haghighi et al. (2010) emphasized the importance of evaluating each model's ability to 
estimate final infiltration rates at different geographic sites. Such location-specific evaluation are 
essential in modeling infiltration to meet the unique requirements of individual regions. Indeed, 
the applicability of infiltration models is not universal, as conducted by Machiwal et al. (2006) 
апа Oku and Aiyelari (2011). Their work affirms that one or more infiltration models may be 
better suited for a particular location, further reinforcing the idea that not all models are 
universally suitable for all soil types and geographic settings. 

The Novel model serves as a valuable tool for generating infiltration data specifically tailored 
to the soils within the Mitidja Plain. The insights gained from this work contribute to the broader 
goal of enhancing water resource management in the Mitidja Plain. The enhanced modeling 
provided by this research offers a more comprehensive understanding of water flow dynamics, 
which is vital for safeguarding the region against flooding and ensuring the efficient use of water 
resources. In summary, the Novel model's capacity to generate infiltration data tailored to the 
Mitidja Plain's soils has far-reaching implications. 


5 Conclusions 


Infiltration models have become one of the most important tools for managing water resource. In 
this work, an attempt was made to validate the predictive accuracy of five infiltration models 
using measurements from a double ring infiltrometer. A comparison of field and predicted 
infiltration rates showed that the Novel model-predicted infiltration rates were much closer to the 
observed data. The NRMSE and KGE results show that the infiltration rate can be predicted by the 
Novel, modified Kostiakov, Horton, Kostiakov, and Philip models. Based on the mean NRMSE 
values, we concluded that the Novel model provided the lowest values, and that the infiltration 
was adequately described by this model. As shown by the results, the Novel model had the lowest 
PBIAS in predicting infiltration rates. Thus, in the absence of measured infiltration data, the 
Novel model can be used to generate infiltration information for the Mitidja Plain and the other 
plains in northern Algeria and North Africa. 


1488 JOURNAL OF ARID LAND 2023 Vol. 15 No. 12 


Conflict of interest 


The authors declare that they have no known competing financial interests or personal relationships that could 
have appeared to influence the work reported in this paper. 


Acknowledgements 


This study was conducted by the Water and Environment Engineering Laboratory (GEE) of the Higher National 
School of Hydraulics (ENSH), and this study is carried out within the framework of the SWATCH project (Prima 
project) funded by the DGRSDT, Algeria. 


Author contributions 


Conceptualization: Amina MAZIGHI, Mohamed MEDDI; Methodology: Amina MAZIGHI, Mohamed MEDDI, 
Giovanni RAVAZZANI, Mouna FEKI; Formal analysis: Amina MAZIGHI, Mohamed MEDDI, Ishak ABDI; 
Writing - original draft preparation: Amina MAZIGHI, Mohamed MEDDI, Hind MEDDI; Writing - review and 
editing: Amina MAZIGHI, Hind MEDDI; Funding acquisition: Mohamed MEDDI, Hind MEDDI; Resources: 
Amina MAZIGHI, Mohamed MEDDI, Giovanni RAVAZZANI, Mouna FEKI, Ishak ABDI, Hind MEDDI; 
Supervision: Mohamed MEDDI, Hind MEDDI. 


References 


Adeniji F A, Umara B G, Dibal J M, et al. 2013. Variation of infiltration rates with soil texture: A laboratory study. International 
Journal of Engineering and Innovative Technology, 3(2): 454—459. 

Angelaki A, Singh Nain S, Singh V, et al. 2021. Estimation of models for cumulative infiltration of soil using machine learning 
methods. ISH Journal of Hydraulic Engineering, 27(2): 162-169. 

Bennett N G, Croke B F W, Guariso G, et al. 2013. Characterising performance of environmental models. Environmental 
Modelling and Software, 40: 1—20. 

Boufekane A, Meddi H, Meddi M. 2020. Delineation of groundwater recharge zones in the Mitidja Plain, North Algeria, using 
multi-criteria analysis. Journal of Hydroinformatics, 22(6): 1468—1484. 

Boufekane A, Yahiaoui S, Meddi H, et al. 2021. Modified DRASTIC index model for groundwater vulnerability mapping using 
geostatistic methods and GIS in the Mitidja Plain Area (Algeria). Environmental Forensics, 23(5—6): 539-556. 

Bouziane O, Mohamed M, Juan Reca C. 2021. Characterization of the water holding capacity of the soils in the Mitidja Plain 
(Algeria) as a basis for the development of the irrigation water use. Journal of Taiwan Water Conservancy, 69(4): 27-39. 

Dahak A, Boutaghane H, Merabtene T. 2022. Parameter estimation and assessment of infiltration models for Madjez Ressoul 
Catchment, Algeria. Water, 14(8): 1185, doi: 10.3390/w14081185. 

Duan R B, Clifford B F, John B. 2010. Field evaluation of infiltration models in lawn soils. Irrigation Science, 29(5): 379-389. 

Esfandiari M, Maheshwari B L. 1997. Application of the optimization method for estimating infiltration characteristics in 
furrow irrigation апа its comparison with other methods. Agricultural Water Management, 34(2): 169-185. 

Farid H U, Mahmood-Khan Z, Ahmad L et al. 2019. Estimation of infiltration models parameters and their comparison to 
simulate the onsite soil infiltration characteristics. International Journal of Agricultural and Biological Engineering, 12(3): 
84—91. 

Golmohammadi G, Prasher $ О, Madani A, et al. 2014. Evaluating three hydrological distributed watershed models: 
MIKE-SHE, APEX, SWAT. Hydrology, 1(1): 20-39. 

Haghighi F, Gorji M, Shorafa M, et al. 2010. Evaluation of some infiltration models and hydraulic parameters. Spanish Journal 
of Agricultural Research, 8(1): 210. 

Iovino M, Angulo-Jaramillo R, Bagarello V, et al. 2017. Thematic issue on soil water infiltration. Journal of Hydrology and 
Hydromechanics, 65(3): 205—208. 

Jain L, Chakma S. 2023. Effects of temperature and slope on the infiltration rate for a landfill surface. Current Science, 124(1): 
94—101. 

Knoben W J M, Freer J, Woods R. 2019. Technical note: Inherent benchmark or not? Comparing Nash-Sutcliffe and 
Kling-Gupta Efficiency Scores. Hydrology and Earth System Sciences, 23(10): 4323—4331. 

Kumar M, Sihag P. 2019. Assessment of infiltration rate of soil using empirical and machine learning-based models. Irrigation 
and Drainage, 68(3): 588—601. 

Lei G H, Fan G L, Zeng W J, et al. 2020. Estimating parameters for the Kostiakov-Lewis infiltration model from soil physical 


Amina MAZIGHI et al.: Estimation and inter-comparison of infiltration models in the... 1489 


properties. Journal of Soils and Sediments, 20(1): 166-180. 

Machiwal D, Jha M K, Mal B C. 2006. Modelling infiltration and quantifying spatial soil variability in a wasteland of 
Kharagpur, India. Biosystems Engineering, 95(4): 569—582. 

Meddi M, Eslamian S. 2021. Uncertainties in rainfall and water resources in Maghreb countries under climate change. In: Leal 
Filho W, Oguge М, Ayal D, et al. Handbook of Climate Change Adaptation. Cham: Springer, 1967—2003. 

Mishra S K, Tyagi J V, Singh V P. 2003. Comparison of infiltration models. Hydrological Processes, 17(13): 2629-2652. 

Ogbe V B, Jayeoba J O, Ode S O. 2011. Comparison of four soil infiltration models on a sandy soil in Lafia, southern guinea 


savanna zone of Nigeria. Production agricultural and Technology, 7: 116-126. 

Oku E, Aiyelari A. 2011. Predictability of Philip and Kostiakov infiltration models under inceptisols in the humid forest zone, 
Nigeria. Kasetsart Journal-Natural Science, 45(4): 604—611. 

Parhi P K, Mishra S, Singh R. 2007. A modification to Kostiakov and Modified Kostiakov infiltration models. Water Resources 
Management, 21(11): 1973-1989. 

Philip J R, De Vries D A. 1957. Moisture movement in porous materials under temperature gradients. Transactions, American 

Geophysical Union, 38(2): 222—232. 

Rasool T, Dar A, Wani M A. 2021. Comparative evaluation of infiltration models under different land covers. Water Resources, 
48(4): 624—634. 

Richards L A. 1931. Capillary conduction of liquids through porous mediums. Physics, 1(5): 318—333. 

Schumann A H. 1998. Infiltration: Introduction. In: Herschy R W, Fairbridge R W. Hydrology and Lakes. Encyclopedia of Earth 
Science. Dordrecht: Springer, 418. 

Sihag P, Tiwari N K, Ranjan S. 2017. Estimation and inter-comparison of infiltration models. Water Science, 31(1): 34-43. 

Singh B, Sihag P, Singh K P. 2018. Comparison of infiltration models in NIT Kurukshetra campus. Applied Water Science, 8(2): 
54. 

Soil Science Equipment | Double Ring Infiltrometer Surface Measure. 2020. SDE Council. [2022-10-22]. 
https://environnement. sdec-france.com/double-ring-infiltrometer-for-soil-science.html. 

Tadj M, Benmouiza K, Cheknane A, et al. 2014. Improving the performance of PV systems by faults detection using GISTEL 
approach. Energy Conversion and Management, 80: 298—304. 

Thomas A D, Ofosu A E, Amuzu E, et al. 2020. Comparison and estimation of four infiltration models. Open Journal of Soil 
Science, 10(2): 45—57. 

Vand A S, Sihag P, Singh B, et al. 2018. Comparative evaluation of infiltration models. KSCE Journal of Civil Engineering, 
22(10): 4173-4184. 

Wang T T, Stewart C E, Ma J B, et al. 2017. Applicability of five models to simulate water infiltration into soil with added 
biochar. Journal of Arid Land, 9(5): 701—711. 

Zakwan M, Muzzammil M, Alam J. 2016. Application of spreadsheet to estimate infiltration parameters. Perspectives in 
Science, 8: 702—704. 

Zeybek M. 2018. Nash-Sutcliffe efficiency approach for quality improvement. Journal of Applied Mathematics and 
Computation, 2(11): 496—503. 


